NEWTONIAN LIMITS OF WARP DRIVE SPACETIMES 
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Abstract. We find a class of warp drive spacetimes possessing Newtonian lim- 
its, which we then determine. The same method is used to compute Newtonian 
limits of the Schwarzschild solution and spatially flat Priedmann-Robertson- 
Walkcr cosmological models. 



Introduction 

Newtonian limits of rclativistic spacetimes are useful tools for understanding 
and interpreting these (see [Ehl91^, _Ehl97^ and references therein). The warp drive 
spacetime, first introduced in |Alc94j . is an interesting example of a relativistic 
spacetime, sheding light on the nature and limits of General Relativity (see for 
instance jLy04 and references therein), but its interpretation as a gravitational 
field is somewhat unclear. It would therefore be desirable to obtain a Newtonian 
limit of this spacetime; this is the problem we address. 

This paper is divided into five sections. In the first section, we quickly review 
the initial value formulation of General Relativity, including formulae which we will 
need in what follows. In the second section, we define Newtonian spacetime and 
its associated potential function, and show that there exist remarkable similarities 
between these spacetimes and their counterparts in Newtonian theory. These sim- 
ilarities are even more impressive in the important case in which the source is a 
dust (with cosmological constant), which we examine in the third section. In the 
fourth section, we show that we can always obtain a Newtonian limit for a New- 
tonian spacetime corresponding to the same potential function; we do this in the 
fifth section for warp drive Newtonian spacetimes, the Schwarzschild solution and 
spatially flat Friedmann-Robertson- Walker cosmological models. 

We shall use the notation and sign conventions in [Wal84| . where latin indices 
represent abstract indices; we will however reserve the latin indices i, j, . . . for nu- 
merical indices referring to the space manifold (i.e., — 1, 2, 3). Whenever 
a tensor field is described by its components, these will always be relative to the 
Eulerian observers' frame defined in the first section. As is usual, we will not worry 
about the vertical position of indices when using Cartesian coordinates. 

1. Synchronized observers 

Let (M, g) be a 4-dimensional Lorentzian manifold, E C Af a spacelike hyper- 
surface and n" the unit normal vector field to this surface. Taking n° as the initial 
condition for a (timelike) geodesic, we can construct a family of free-fall observers 
(which we shall call Eulerian observers). The proper time of each observer along 
the geodesic (starting from S) defines a time function t in an open neighbourhood 
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of S. Shrinking E if necessary, we can assume that t fohates some open set U by 
spacehke surfaces Y,tg = i^^(to), with Sq = S, diffeomorphic by the flow of n". 
As is well known, the hypersurfaces Stp are orthogonal to n", and hence integrate 
the distribution of hyperplanes defined by the local notions of simultaneity of the 
Eulerian observers; we then say that the Eulerian observers are synchronized - see 
[LL971 V 

The integral hues of determine a natural projection tt : [/ ^ E. Hence 
local coordinates on an open subset V C 'S yield local coordinates (t,^') on 
TT^^{V) C U. In these coordinates, 

g = -dt (g) dt + (t) d^* (g) d^^ . 

Notice that 

can be thought of as a (time-dependent) Riemannian metric on S. We shall call 
(E,7(t)) the space manifold at time t. 

It is a simple matter to compute the components of the Riemann curvature 
tensor R^^^,^ and of the Ricci tensor Rab of (M, g) on these coordinates: they are 

(1) R,,,^ ^~^K\-KuK'=- 



and 



Rm = -^^K\-K,,K'^- 
RQr^-D,K\ + DjK'- 



where R^ji ™ = R^ji and Rij — Rij[t) are the components of the Riemann and 
Ricci tensors of the Levi-Civita connection D = D{t) of 7(i), and 

In these calculations, one must bear in mind that raising and lowering indices with 
the space manifold metric ^it) does not commute with since 

d d 
— 7„- = 2X„-, — 7'J = -2K'^. 
dt ' ^ ^' 9t ' 

For instance, equation ([T]) can be rewritten as 

(2) = -| h''Ku) - KuK'' ^ + KuK^'- 

Again, 

X(t) = K,j{t)dC®d^^ 

can be thought of as a (time-dependent) symmetric tensor on S. On the other 
hand, both 7 and K can be thought of as tensor fields on U, by considering j{t), 
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K{t) to be tensor fields on each submanifold St C J7 and then extending them 
trivially to the direction given by n": 

Jab = gab + naHb] 

Kab = -j^£nlab = ■^£ngab = ^ (^a"'' + ^fc"a) • 

Notice that Kab yields the extrinsic curvature of the family Sj. 

2. Newtonian spacetimes 

If all the hypersurfaces are flat, then one can easily find coordinates {t, x') in 
U such that (x' ) are Cartesian coordinatefQ on : 

(3) g^~dt(g>dt + 5ij {dx' ~ v'dt) <S> {dx^ - v^dt). 

Notice that on these coordinates one has 

^a ^ ^ 

dt dx'- ' 

therefore, v = f^^r can be intrepreted as the 3- velocity of the Eulerian observers 
on a fixed Euclidean background. Since the coordinates are are only defined up 
to a (time-dependent) Euclidean isometry, v is only defined up to the addition of 
an arbitrary (time-dependent) Euclidean Killing vector field u}{t) x x + xo(t), which 
is, of course, irrelevant when computing 

K — 2'^n7 ~ + DjVi)dx^ (g) dx\ 

If V can be chosen to be a gradient^, v = V?/', then the Eulerian observers will 
follow Newtonian free-falling motions on the potential 

as in this case DiV^ = Djv'^ and hence 



(4) - D,cj) = - (A^) + vW^v^ = ^'D,v' ^ 7^ + (v • V)v = -V0. 

ot at ot 

Definition 2.1. We say that a metric of the form ([3]) with v — Vifj defines a 
Newtonian spacetime with potential <P — ~ 5^^- 

Proposition 2.2. The components of the Riemann and Ricci tensor for a Newto- 
nian spacetime are 

(5) Ro^o' = DW,q^- 

(6) i?,^o' = 0; 

(7) R^^r = KuK"^, - K,iK^,- 



^ Using for instance the exponential map on these hypersurfaces. The coordinates are some- 
times called Eulerian coordinates, whereas ^* are known as the Lagrangian coordinates. It would 
be perhaps more reasonable to call the Eulerian observers 'Lagrangian observers', as it is their La- 
grangian coordinates which remain constant; nevertheless, we shall keep with the common General 
Relativity usage. 

^Note that this will fix v up to a (time-dependent) additive constant xo(t). 
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Roo = 



Ro^ = 0; 

(8) R^J = -D.Djcj) - KaK'j + K\K,, . 

Proof. For a Newtonian spacctime Kij = Div^ . Equation fl]) then yields 



dt 

from which 

d_ 

dt 
d 



(9) {K\) = + = - - if'^ . 

Using ([2]), the result follows. □ 



Newtonian spacetimes share a number of common features with their Newtonian 
counterparts: 

Proposition 2.3. The potential (f> satisfies the Poisson equation 

V20 = 47r(Too+7'''7^u)-A. 
here Tab is the energy-momentum tensor and A is the cosmological constant. 
Proof. This is just the (00) component of Einstein's equation 

(10) Gab + Agafc = 87rT„fc Rab = StT [Tab - ^T"^ 9ab^ + Agab- 

□ 

Proposition 2.4. Free-falling particles on a Newtonian spacetime follow, to first 
order on their velocities with respect to the Eulerian observers, Newtonian free- 
falling motions with respect to the potential 4>. 

Proof. Using t as the parameter, the Lagrangian for free-falling motions on a New- 
tonian spacetime is 



To second order on the velocity with respect to the Eulerian observers this is 

2 ^ \ dt J \ dt 

The Lagrangean for a Newtonian free-falling motion with respect to the potential 
is 

, _ 1 dx^ dx^ 

~ 2 "^~dt~dt ~ ^' 
Using the definition of 0, we see that 

and hence L and L' determine the same motion equations. □ 
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Proposition 2.5. // a Newtonian spacetime is flat then the potential is an affine 
function. 

Proof. This is immediate from equation ([3]). □ 

Surprisingly, the converse is not true. For instance, taking v = jx one obtains 
= and K = jl, and hence 

(11) (i?..)=(2 4, 

This Newtonian spacetime corresponds to a cosmological solution with a negative 
pressure perfect fluid which nevertheless satisfies (just) both the strong and the 
weak energy conditions. 

3. Dust Newtonian spacetimes 

To get a closer analogue to the Newtonian situation, we consider the case when 
the source of the gravitational field described by a Newtonian spacetime is a dust 
comoving with the Eulerian observers. 

Tab = pnaUb 

(this includes vacuum solutions; see also [Tra66| ) . In this case, the Poisson equation 
in Proposition 12 . 31 is simply 

(12) V20 = 47rp-A, 

which is exactly the Newtonian Poisson equation with mass density p and cosmo- 
logical constant A. We also have 

Proposition 3.1. For Newtonian dust fields with rest density p, the continuity 
equation 

|+V.(pv)=0 

holds. 

Proof. Since n° is geodesic, one has 

(13) V'Tafc = 0^n^VbP + p^bn^ = 0^ £nP + ptrK = 0, 
which can be recast as the usual continuity equation: 

|j + (v-V)p + pV-v = 0^^+V - (pw) = 0. 
at ot 



□ 



Dust fields Newtonian spacetimes behave satisfactorily in another respect: 



Proposition 3.2. If the potential associated to a Newtonian dust field is an affine 
function then the solution is flat. 

Proof. Suppose that 770 = 0. Then = 0, and hence from Poisson's equation 
(dH) one has Anp = A. 

If A ^ 0, one has from (O that tvK = 0. Equation ^ then yields tr = => 
K = 0, and it is clear from Q and (O that the solution is flat. 

If A = 0, one has p = and Einstein's equation pil| yields i?^^ = 0. It follows 
from (O that K'^-{ti K)K = 0. From we see that R^^^ ' = 0; since ii^, ™ has the 
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same symmetries of the Riemann tensor of a 3-dimensional Riemannian manifold, 
we conclude that i?^, " = 0. Finally, ^ and ^ yield Eg-g^ = E.^o' = 0. □ 

4. Newtonian limits of Newtonian spacetimes 
To construct a Newtonian limit of a Newtonian spacetime we consider the family 

(14) g{\) = ~^dt ®dt + Sij{dx' - v'dt) ® {dx^ - v'dt) 

A 

(see |Ehl91j . [EhlQTj V 

Proposition 4.1. Family ^4\ l a Newtonian limit corresponding to a gravita- 
tional potential (p. (In other words, Newtonian spacetimes always have Newtonian 
limits with the same potentials). 

Proof. This is immediate from Proposition [5131 which shows that the leading terms 
of the Christoffel symbols for the Levi-Civita connection of a Newtonian spacetime 
are the Christoffel symbols for the Cartan connection corresponding to a gravi- 
tonal potential (p. As confirmation, we compute the components of the Riemann 
curvature tensor for this family: 

As A ^ these clearly converge pointwise to the components of the Riemann tensor 
of the Cartan connection corresponding to a gravitational potential (j). □ 

Note that in general both the Newtonian spacetimes and their Newtonian 
limits are unphysical, in the sense that they do not come with a corresponding 
matter model generating the field (this must be so to accomodate the warp drive) . 
In the particular case of a dust, however, they do yield a genuine physical Newtonian 
limit in the sense of |Ehl91| . |Ehl97| . 

5. Examples 

5.1. Newtonian warp drive. A warp drive spacetime containing a warp bubble 
moving with speed u{t) is given by a metric of the form ^ with 

{0 if r < ro — £ 

-u{t)-£r ifr>ro + e 

where r^ = (a;^)^ + (x^)" + (x^)^ and tq > e > (this is in the frame where the 
bubble is at rest; see |Nat02| for deatils). If v = Vi/;, this will also be a Newtonian 
spacetime. Choosing 

ij = -u{t)f{r)x\ 
where / : M ^ [0, 1] is a smooth function satisfying 

/(r) = if r < tq — e 
/(r) = 1 if r > ro + e 

it is easily seen that 

d 

V = Vtb ~ —uf-——^ — uf'x^\7r 
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satisfies the above requirements. The corresponding potential is 



2 ^ 



dt 2 



v2=^./xl-^ + +2//' 



This potential contains two terms: an acceleration term, which is only present 
when the velocity u{t) is changing, and a steady term, which must be present 
whenever the warp drive is on {u{t) ^ 0). 

One has = for r < tq — s and (j) = iix^ ~ for r > tq + e. Thus outside the 
bubble wall region one has vacuum and, if the bubble is not accelerating, vanishing 
gravitational fields. An uniform gravitational field of magnitude u is responsible 
for accelerating the bubble (thus guaranteeing zero spacetime curvature outside the 
bubble). 

By ordinary potential theory, one knows that cannot be identically zero, 
and that moreover 

l|x-xo|| 

for all xq satisfying ||xo|| > rg+e. Therefore, there will be points on the warp bubble 
wall where V^(/) < 0. At these points, one has Rqq < 0, and the strong energy 
condition is violated (as it we know it must be - see [Nat02j ): in the corresponding 
Newtonian limit, this translates as the mass density generating the gravitational 
field being negative. 

5.2. Schwarzschild solution. It is well-known that the Schwarzschild metric can 
be put in the so-called Painleve- Gullstrand form 

g = -dt ®dt + Sij dx* - 1 dt \ 1^ \dx^ ~ 1 dt 

y r2 y y r2 J 

(see |Vis04) ) . Since the vector field v with components 

i _ {2M)ix' 
^ I 

7" 2 

is clearly a gradient, we immediately see that this corresponds to a Newtonian 
spacetime with potential 

which is the Newtonian potential of a point mass M, as one would expect. 

Notice that the horizon corresponds to |jv|| = 1; this can be seen as a consequence 
of the fact that the condition for a timelike curve (using the time coordinate t as 
the parameter) is ||^ ~ v|| < 1, and explains why the Newtonian formula for the 
Schwarzschild radius coincides with its General Relativity counterpart. 

5.3. Spatially fiat Priedmann-Robertson- Walker models. Another simple 
but instructive example is to consider the Newtonian cosmological model gener- 
ated by the choice 

a = u{t)x' ^ D ' ^ ' 

corresponding to the potential 

d (u{t)r^\ 1 



dt\ 2 2 



VjV — — 7i + 7i 



2^ 
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We have Hcf) = — (u + u )I and K = ul: consequently, 

f-i3u + 3u^) \ 

^^t^'^'-y {u + 3u^)l) ^ ^ \ -{2u + 3u^)I 

(see also pTjl ). Einstein's equation PH)) for a dust then yields 

JSu^ - A = 87rp 
[2u + 3u2 - A = 

One can solve the second equation for u{t) and then obtain p{t) from the first. 
These have the following nontrivial physical solutions (the last three with time 
coordinate chosen so that the big bang occurs at t = 0): 

(1) De Sitter universe (A > 0): 

A , Ar2 

3' ^ = °' '^ = -3y 

(2) Spatially flat Priedmann-Robertson- Walker model with A = 0: 

2 1 , 2 r2 



(3) Spatially flat Priedmann-Robertson- Walker model v^^ith A > 0: 

A . /I rr^\ A (a 



a/— cotanh I -A/3At I , p 
V 3 V2 / 



Stt sinh^ ( i VSAt) \ 6 sinh^ ( i VSAt 



(4) Spatially flat Priedmann-Robertson- Walker model v^rith A < 0: 



Stt sin^ 



(i^sRi)' \6sin2(iV3Ri) 



These are good examples of how Newtonian limits (described by the time- 
dependent potentials above) can give insight int the nature of the solutions. See 
[NST99| for more examples of Newtonian cosmological solutions of Einstein's field 
equations. 
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